Mobility of thorium ions in liquid xenon 



in 
o 
o 



(N 



5h 

43 



K. Wamba, 1 C. Hall, 1 M. Breidenbach, 1 R. Conley, 1 A. Odian, 1 C.Y. Prescott, 1 P.C. Rowson, 1 
J. Sevilla, 1 K. Skarpaas, 1 R. DeVoe, 2 Z. Djurcic, 4t W.M. Fairbank Jr., 5 G. Gratta, 2 M. Green, 2 K. Hall, 6 
M. Hauger, 6 S-C. Jeng, 5 T. Koffas, 2 A F. LePort, 2 D. Leonard, 4 J-M. Martin, 6 R. Neilson, 2 L. Ounalli, 6 
A. Piepke, 4 D. Schenker, 6 V. Stekhanov, 3 J-L. Vuilleumier, 6 S. Waldman, 2 P. Weber, 6 J. Wodin 2 
1 Stanford Linear Accelerator Center, Menlo Park CA 94025, USA 

2 Physics Department, Stanford University, Stanford CA 94305, USA 

3 Institute for Theoretical and Experimental Physics, Moscow, Russia 

4 Department of Physics and Astronomy, University of Alabama, Tuscaloosa AL 35487, USA 
5 Physics Department, Colorado State University, Fort Collins CO 80523, USA 
6 Institut de Physique, Universite de Neuchatel, Neuchatel, Switzerland 

(Dated: February 2, 2008) 

We present a measurement of the 226 Th ion mobility in LXe at 163.0 K and 0.9 bar. The result 
obtained, 0.240±0.011 (stat) ±0.011 (syst) cm 2 /(kV-s), is compared with a popular model of ion 
transport. 
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I. INTRODUCTION 

Free ions in the liquefied noble gases at low concentra- 
tions represent a theoretically straightforward condensed 
matter system. This is largely because of the liquid's 
atomic nature, implying relatively few microscopic inter- 
nal degrees of freedom. 1 Despite this fact, the experi- 
mental study of these systems has long been confounded 
by the difficulty in reproducibly generating and trans- 
porting ions of known identity while controlling for the 
effects of impurity ions. Accordingly, there have been 
relatively few transport studies for ions in the liquefied 
noble gases other than LHe reported in the literature to 
date. Furthermore, studies that deal only with molecular 
ions ) 1 ' 2 ! 3 ! 4 ! 5 are more common than those that deal with 
monatomic speciesAi There is therefore a clear need for 
additional data on monatomic ion transport particularly 
in the heavy noble gases. In the present work we consider 
the system of monatomic Th ions in LXe. 

This research also has important applications to the 
physics of ionizing particle detectors, especially in the 
area of rare event searches that use a liquefied noble 
gas as the detection medium&^iiS For these systems, the 
transport properties of Th and U ions are of concern be- 
cause these elements in trace amounts represent a large 
source of background radioactivity. It is therefore de- 
sirable to understand how electric fields might be used 
to draw these ions out of the detection medium^! For 
LXe-based detectors, Cs and Ba ions are also interest- 
ing because their respective separation from the detec- 
tion medium and subsequent identification can greatly 
enhance the sensitivity to solar neutrinos and to double- 
beta decay£*i The motivation for the present study is 
work involving Ba ion retrieval from LXe as part of R&D 
for a large double-beta decay experiment^ 

In this paper, we report on the measurement of the 
mobility of 226 Th ions in liquid Xe at 163.0 K and 0.9 
bar pressure. We also compare our results and those of 



other workers with the theoretical calculations of Hilt, 
Schmidt and Khrapak. 2 



II. EXPERIMENT 

Our system, schematically shown in Fig. ^ consists of 
a vacuum-insulated cell which we fill with about 10 cm 3 of 
LXe. The cell is cooled by a cold finger that forms a ther- 
mal link with a liquid nitrogen reservoir. The tempera- 
ture is maintained at 163.0±0.1 K by a resistive heater. 
The cell, built almost entirely out of UHV-compatible 
materials, is baked at 400 K for several days, after which 
it is pumped to a 7.0xl0~ 8 mbar vacuum by means 
of two turbomolecular pumps. When filling the cell, we 
pass research-grade Xe (which, as supplied, is quoted as 
99.999% pure) through a hot Zr getter^ which is rated 
to remove reactive impurities to the <1 ppb level. These 
precautions are taken to minimize any effect of impuri- 
ties that might influence ion transport characteristics. 

Th ions are produced by a 230 U source mounted 
inside the cell. The 230 U was obtained as a de- 
cay product of 230 Pa, in turn produced in the re- 
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>g 30 Pa+3n on natural thorium at a 
cyclotron^ After irradiation, the 230 Pa was left to (3- 
decay for a few of its 17.4 day half lives, allowing some 
build-up of 230 U, which was then chemically separated 
and electroplated onto a platinum disk. Fig. [2Jshows the 
decay scheme for the final 230 U source. 

226 Th ions are formed as a result of the 20.8 day a 
decay of 230 U. The average recoil energy of ^100 keV ex- 
tracts 226 Th atoms from the surface in a highly ionized 
state. As the 226 Th ion enters the liquid it exchanges 
charge with the surrounding Xe atoms and reaches equi- 
librium in some final charge state. Although the energet- 
ics of the charge exchange in principle favor a final state 
containing Xe + and Th + (Ref. 0), a large number of the 
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FIG. 1: Schematic of the probe cell. 



emerging Th ions can also recombine with free ionization 
electrons that arise from the original a decay£ This effect 
presents a major problem as it greatly suppresses the ion 
yield. 

We collect ions on a stainless steel probe which we 
lower into the liquid. To collect ions we place the probe 
at ground potential while applying a positive bias to the 
source electrode. The tip of the probe is fashioned into 
a 5 mm x 16 mm, 2 mm thick rectangular stainless steel 
plate, ensuring a spatially uniform electric field in the 
vicinity of the source. The uniformity of the field is veri- 
fied by constructing a computer model of the electrostat- 
ics with the MAXWELL software packaged 

The ions that land on the probe are detected by re- 
tracting the probe vertically to a counting station, lo- 
cated above the liquid surface. It consists of a 50 mm 2 
active area ion-implanted Si detector— that is read out 
using standard pre- and post-amplification electronics. 
The ion yield is determined by counting the a particles 
emitted in the 30.6 min 226 Th decay, along with those 
produced by the subsequent a decays of 222 Ra, 218 Rn, 
and 214 Po (see Fig. EJ. 

It should be pointed out that in addition to the 226 Th 
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FIG. 2: Decay scheme for 230 U. Decays with branching ratios 
less than 5% are not shown. 



ions, the 230 U source also produces ions arising from the 
other a decays in the chain. These can also drift to the 
probe and be detected at the counting station. How- 
ever, because of their considerably shorter lifetimes, they 
present a negligible background to the 226 Th signal pro- 
vided some short time elapses between the probe extrac- 
tion and the counting. In addition, the low angular ac- 
ceptance of our a particle detection system (~5%) en- 
sures that the probability of detecting more than one a 
particle from the 4-a decay sequence of a given Th atom 
is negligibly small. 

We measure the ion mobility in liquid Xe by deter- 
mining the ion transit time from source to probe at a 
given probe distance and electric field. This is achieved 
by biasing the source with an alternating potential pro- 
duced by a computer driving a high voltage amplifier. 16 
The system is set to produce a square wave switching be- 
tween Vdrift and V sup with respect to the probe, where 
Vdrift can be varied from 1.0 kV to 4.0 kV, and V sup is 
fixed at -4.0 kV. Vdrift is set for a time tdrift ~ 1 s, very 
long compared to the ^20 fxs switching time of the high 
voltage system. During each full cycle, the voltage is re- 
versed for a time t sup . When this occurs, the ions that do 
not reach the probe during tdri ft are drawn back towards 
the source. This occurs when tdrift is too short for the 
ions to complete the journey to the probe. We maximize 
our total ion integration time by having t sup < tdri ft and 
\Vdrift\ < \V sup \A During each data taking run, Vdrift is 
set to a particular value and the field is switched for 1800 
s. The constant duration of collection runs ensures that 
the amount of ion loss due to radioactive decay during 
collection is uniform. We then measure the ion yield at 
the counting station over a period of 2 hours to allow suf- 
ficient time for the collected 226 Th ions to decay before 
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the next run. Data are taken at four different electric 
fields, using four source-probe voltages and two source- 
probe separation distances. 

We determine the background by performing 2-hour 
counts between data runs in which no ions are collected. 
In these runs the source-probe voltage remains reverse- 
biased. The observed rate during the background runs is 
averaged over all runs, giving the mean background rate 
and an associated statistical error. The resulting value of 
3.16±0.73 counts per 2-hour run is subtracted from the 
number of counts observed at each tdri ft ■ When ions are 
collected, the typical rate is —30 counts per run. Sources 
of background include cosmic rays and neutral Th atoms 
that evaporate from the surface of the liquid and diffuse 
through the gas. 

For the purpose of characterizing the system a test in 
Xe gas was also done. The cell was filled with 0.2 bar 
of room-temperature Xe gas and the probe was alter- 
nately brought close to the 230 U source and raised to the 
counting station, where ions were respectively collected 
and counted. This was repeated many times over several 
hours to build up counting statistics. The electric field 
between source and probe for ion collection was roughly 
-200 V/cm for this test. Shown in Fig. are the a 
spectrum and time histogram for the collected Th and 
its decay daughters. For comparison, the 230 TJ spectrum, 
as measured in vacuum at LLNLjli is also shown. An 
important conclusion from the gas test was that ion re- 
combination with free ionization electrons is a factor of 
—200 times more efficient in LXe than in gaseous Xe. 




Time (s) 

FIG. 3: (a) a spectrum of the 230 U source, taken under 
vacuum!— ■ The 230 U peaks are apparent at 5.82 MeV and 5.89 
MeV, while 226 Th is identified by the peaks at 6.23 MeV and 
6.34 MeV. 222 Ra, 218 Rn, and 214 Po correspond to the peaks at 
6.56 MeV, 7.13 MeV, and 7.69 MeV respectively, (b) a-decay 
spectrum for ions that are collected on the probe in 0.2 bar of 
ambient Xe gas. Here, the peaks are broadened and shifted 
to lower energy as a result of losses in the Xe gas. (c) a- 
decay time histogram for the sum over all counts in the alpha 
spectrum shown in (b). The best fit is the superposition of two 
exponentials with half-lives of 1826±106.7 s and 38.2±2.97 s. 
These match the 226 Th and 222 Ra half-lives, whose published 
values are 1830 s and 38 s respectivelyi- 



III. RESULTS 

Shown in Fig. 0|are plots of the ion yield vs t&rift- For 
each data point, the ion yield is determined by dividing 
the number of detected a particle counts by the num- 
ber of cycles, giving the ion yield per cycle. The average 
background in these units is <0.002, and the statistical 
error bounds on each data point were computed using 
the Feldman-Cousins prescription^ We find that the ion 
yield is consistent with zero (after background subtrac- 
tion) when tdrift is short, while at longer times, the ion 
yield per cycle increases linearly with tdri ft ■ The abscissa 
of the join point between these two behaviors is the ion 
transit time. We fit the ion yield data to the following 
two-parameter piecewise linear function, 

Y = tdrift < ttr 

Y = m(tdrift — ttr) tdrift > t tr 

where Y is the ion yield per cycle, m is the slope (ion 
yield per cycle per unit tdrift), and t tr is the ion tran- 
sit time. The value of t tr yielded by the fit is used to 
calculate the ion drift velocity. To perform the fits, and 
to calculate the statistical error on tt r in each case, we 
use an algorithm due to D.J. HudsonSi that is designed 
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FIG. 4: Ion yield Y as a function of drift time tdrift for the four electric fields 2.41 kV/cm, 3.93 kV/cm, 4.86 kV/cm, and 6.42 
kV/cm. Error bars shown are statistical only. We take the transit time in each case to be the knee point in the fit. 



TABLE I: Electric fields (E), voltages (V drif t), and drift dis- 
tances (d) used in this experiment, showing their associated 
systematic errors. The measured transit times (tt r ) and ve- 
locities (u) are the values returned by the fit. For these last 
two quantities, only statistical errors are given. 



d 

(xlO cm) 


Vdrift 

(xlO kV) 


E 

(kV/cm) 


ttr 

(xlO s) 


V 

(x 10 cm/s) 


3.87±0.08 
5.87±0.08 
3.87±0.08 
5.87±0.08 


9.32±0.14 
23.10±0.24 
18.88±0.10 
37.71±0.32 


2.41±0.06 
3.93±0.07 
4.86±0.11 
6.42±0.11 


6.96±gJl 
6 64+ 1 01 

o 9n +0.26 
°- zu -0.25 
o o n +0.40 
°- ou -0.35 


c kc+0.27 
o.ou_ 0A6 

19 1+0.98 
lz - 1 -0.94 



specifically to handle the difficulties that arise when per- 
forming a piecewise fit. One peculiarity of such a fit is 
that the error bounds are almost always asymmetric, — 
The results of the fits are summarized in Table [I] 

Fig. shows a plot of the computed ion drift veloc- 
ity as a function of electric field. The data are consis- 
tent with the linear relation v = /j,E, where v is the ion 
drift speed, /i is the mobility, and E is the electric field. 
From a fit to the data in Fig. [3] we find an ion mobil- 



ity [i =0.240±0.011 cm 2 /(kV-s). This fit is obtained by 
forcing the straight line to go through the origin and it 
yields x 2 /ndf=2.02/3.A fit not constrained to go through 
the origin results in \i = 0.268 ± 0.031 cm 2 /(kV-s), a 
y-intercept of -0.095±0.088 cm/s, and x 2 /ndf=1.14/3. 
This is fully consistent with the original result. 

The mobility and its statistical error are determined 
largely by the data taken at the lowest electric field. This 
point is more precise for two reasons. First, the ion tran- 
sit time at this electric field is the longest of the four 
datasets, meaning that the fractional error on the ve- 
locity is the smallest. Secondly, in this dataset the ion 
transit time happens to be within ^0.01 seconds of the 
nearest tdri ft at which we collected data, and this results 
in a precise measurement. While, in principle, better ac- 
curacy could be obtained by repeating the measurements 
and always taking a point near this optimal time, the 20.8 
day half life of the 230 U source makes it impractical to 
repeat the entire experiment. It is also true that the data 
from which the first point is determined was recorded ear- 
lier in the lifetime of the 230 U source, resulting in higher- 
statistics. 

Systematic uncertainties on the mobility derive from 
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H = 0.240± 0.011 kV/cm-s 



TABLE II: Ion mobilities in LXe, comparing a model due to 
Hilt, Schmidt and Khrapak (HSK) 2 to the experimental data. 
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FIG. 5: Global fit to the ion drift velocity as a function of 
electric field for the entire data set. 



small mechanical vibrations of the probe when the elec- 
tric field is switched, and from the accuracy of the voltage 
gain in the amplifier. These contributions are estimated 
to give a combined error of ± 0.011 cm 2 /(kV-s). Table 
Q] shows the absolute errors on the voltages, drift dis- 
tances and electric fields. The systematic error on the 
ion mobility was computed by shifting the electric field 
and drift distance values by their corresponding errors 
and re-doing the fit. 



IV. DISCUSSION 

To help us understand our result, we have considered 
a transport model due to Hilt, Schmidt and Khrapak 2 
(henceforth referred to as HSK). In this model the mi- 
croscopic properties of the ion and its detailed coupling 
to the liquid atoms are neglected, and the ionic mobility 
is determined solely by the equation of state of the liquid 
in question. The ion is modeled as a point charge whose 
electric field exerts a pressure on the liquid of sufficient 
magnitude for a microscopic ice layer to form around the 
ion. The ion and its attached cluster of Xe atoms then 
behave as a positively charged ball of ice whose radius is 
determined by the distance at which the electric field is 
no longer of a magnitude sufficient to compress the liq- 
uid into the solid stated. In this theory, no distinction 
is made between the mobilities of ions of different iden- 
tities, and the drag suffered by the ion and its attached 
cluster is computed from the bulk viscosity of the liquid. 

Shown in Table ITT1 are some selected experimental mo- 
bilities for atomic ions in LXe. From these results it is 
clear that the HSK model does quite well in predicting 
the overall scale of the observed mobilities. However, 
the significant difference between the 208 T1 + and 226 Th + 
numbers suggests that these species differ enough micro- 
scopically that the HSK model, when applied to them, 



Ion 


T (K) 


P (bar) 


Mobility (cm 2 kV" 1 s" 1 ) 
HSK Expt. 


208rpj+ 

226 Th + 
Ba+ 


163.0 
163.0 
168.0 


0.9 
0.9 
1.2 


0.31 
0.31 
0.33 


0.133±0.004 a 
O^OiO.OlliO.Oll 6 

n 91 -,+0.020 c 
u - z±1 -0.012 



°Ref. |3 
'This work. 
c Ref. 



begins to break down. 

It should be pointed out that all of the calculations 
assume an ionic charge state of +1, which may not be 
correct. This is based on the argument that the energet- 
ics of charge exchange favor such a final stated but this 
assertion has not yet been convincingly verified experi- 
mentally. Nonetheless, we find that assuming an ionic 
charge other than +1 does not improve the compatibil- 
ity between the theory and the data for any of the cases 
considered here. 

Overall, it can be concluded that HSK provides a rea- 
sonable rough approximation of the ion mobility for the 
systems under consideration. However, it remains to be 
seen whether the model produces the correct temperature 
dependence in all cases. In addition, it is worth under- 
standing the situations in which the ion identity becomes 
important because this would offer clues as to how HSK 
could be simply improved. 

A theory put forward by Davis, Rice and Meyer (hence- 
forth referred to as DRM) might offer some clues as to 
how this might be accomplished. 1 These authors have de- 
rived an expression for the positive ion mobility in which 
the internal structure of the ion and its interaction with 
the liquid are taken into account. In this model the drag 
force experienced by the ion is determined by the micro- 
scopic interaction potentials rather than the bulk proper- 
ties of the liquid. It is conceivable, at least in principle, 
that as the ion-liquid interaction increases in strength, 
bound states between the ion and several liquid atoms 
can occur, particularly as the temperature of the liquid is 
lowered. Nevertheless, the treatment reported in Ref. Q] 
is specific to conditions in which this effect is not present, 
i.e., the ion- liquid interaction is too weak and the liquid 
temperature is too high. Indeed, the DRM model has 
only been compared with experimental data from sys- 
tems involving molecular ions in liquids at temperatures 
well above the liquid triple point, namely He2 + in LHe 
Ref. H, Ar 2 + in LAr (Ref. 0), and Kr 2 + in LKr (Ref. 

While these systems are satisfactorily described by 
the model, it is clear that some modification will be nec- 
essary in order to treat the case of atomic ion mobilities 
close to the triple point. This is because the atomic ions 
under consideration have interaction strengths that are 
tens of times strongesSi^ than those of the molecular 
ions previously considered. 
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To investigate these details experimentally, one could 
perform an accurate measurement of the temperature de- 
pendence of the ion mobility in several different ion-liquid 
systems. This would require multiple, well-understood 
and well-controlled sources of ions of known identity. A 
complementary approach would be to use the 226 Th ion 
source in several different liquefied noble gases at various 
temperatures. 



when used to calculate ion mobilities, however it cannot 
account for the significant difference between our result 
and the one measured for 208 T1+ by other workers. A 
more detailed study of this effect may help to suggest 
how the HSK model can be modified to account for 
different ion identities. A possible route might be to 
generalize the DRM treatment to the case of atomic ions 
that form bound states with liquid atoms. 



V. CONCLUSION 

We have measured the 226 Th ion mobility in LXe at 
163. OK and standard pressure and have found it to have 
a value of 0.240±0.0f f (stat) ±0.011 (syst) cm 2 /(kV-s). 
The HSK model gives the correct order of magnitude 
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